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DIMENSION REDUCTION FOR ROTATING BOSE-EINSTEIN 
CONDENSATES WITH ANISOTROPIC CONFINEMENT 

FLORIAN MEHATS AND CHRISTOF SPARSER 


Abstract. We consider the three-dimensional time-dependent Gross-Pitaevskii 
equation arising in the description of rotating Bose-Einstein condensates and 
study the corresponding scaling limit of strongly anisotropic confinement poten¬ 
tials. The resulting effective equations in one or two spatial dimensions, respec¬ 
tively, are rigorously obtained as special cases of an averaged three dimensional 
limit model. In the particular case where the rotation a:xis is not parallel to the 
strongly confining direction the resulting limiting model(s) include a negative, 
and thus, purely repulsive quadratic potential, which is not present in the origi¬ 
nal equation and which can be seen as an effective centrifugal force counteracting 
the confinement. 


1. Introduction and main result 

We are interested in the dimension reduction problem arising in the descrip¬ 
tion of rotating Bose-Einstein condensates with strongly anisotropic confinement 
potential. In physics experiments such potentials are used to obtain effective one¬ 
dimensional (called cigar-shaped) or two-dimensional (called pancake-shaped) con¬ 
densates which, among other features, exhibit different stability and instability prop¬ 
erties than the usual three dimensional case (for a general introduction to the physics 
of Bose-Einstein condensates, see, e.g., [16, 17]). The present work aims to give a 
rigorous justification to the use of these approximate lower-dimensional models. In 
comparison with earlier studies in the mathematics literature, see [2, 3, 4, 5, 6], 
the main novelty in our work is the presence of an additional angular momentum 
rotation term, whose strong interaction with the confinement will, in general, result 
in a nontrivial effect within the limiting model obtained. 

The starting point of our investigation is the three-dimensional Gross-Pitaevskii 
equation, describing the Bose-Einstein condensate in a mean-held approximation, 
cf. [14, 15, 17]. Rescaled into dimensionless form (see, e.g., [6]) and in a rotating 
reference frame (which is customary used throughout the literature), this nonlinear 
Schrodinger equation (NLS) reads 

idtfi = fi + in■{■K^V)fi + (1.1) 
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with an initial data = 0,x) = and where Q = 

is a given vector of describing the rotation axis. Here, the exponent of the 
nonlinearity is assumed to be 1 < u < 2. In particular, this means that the 
nonlinearity is i?^(M^)-subcritical, cf. [9]. Note that this includes the cubic case 
(7 = 1, which is the physically most relevant situation. The space variable x G 
splits into X = {x,z) G x M, with x = (a;i,X 2 ). Finally, we assume that e G (0,1] 
is a small parameter describing the anisotropy within the confining potential. In the 
following we shall be interested in the limit e —>■ 0 for solutions to (1.1). Note that we 
thereby assume that the initial wave function is already confined at the scale epsilon 
in the z-direction (an assumption which is consistent with the asymptotic limiting 
regime considered) and such that its total mass \\'>p{t = 0, Ollia = Ij uniformly in s. 
Let us rescale the variables. We set 

x'= X, z' = -, 'tp‘^{t,x',z') = (t,x',ez') 

and assume that (3^ = Ae®^, where A G M is fixed. We are thus in a weak interaction 
regime similar to [5, 6]. Under this rescaling the NLS becomes (dropping the primes 
in the variables for simplicity) 

- ^1x2) dzip^ + - ^zLzi^^ 

- i£z {9.idx2 - n25xi) (1.2) 

with i/j^^t = 0,x,z) = ipo{x,z). Here, and in the following, we denote 
1 ^2 |^|2 

T~lz — 2 ^^ 2 ’ — ^^2^X1 ixidx2 ! 

where is the angular momentum operator associated to a rotation around the 
negative z-axis. 

In order to get rid of the singular rotation term proportional to £~^ in (1.2), we 
shall invoke the following (unitary) change of unknown. Setting 

X, Z) = e^ez(niX2-Q2Xi).^e^^^ 3 ) 

we obtain 

idtu^ = ^TizU^ + T-ixU^ — - {^2X1 — ^1X2)^ — VLzLzU^ + 

3^2 " (1-4) 

H —— (flf + nl) ~ (flixi + ^ 2 X 2 ) zu^ + 2i£z{^2dxx — ^idx2)'u‘^ 1 
subject to initial data 

u^{t = 0 ) = u^ = (^ 5 ^ 

Note that in (1.4), the only singular term left is £~^'Hz- We consequently expect that 
by filtering the associated rapid oscillations the function u^{t) will converge 

to some finite limit ^{t), as e —)• 0. 

A suitable functional framework for the analysis of our problem is the scale of 
Sobolev spaces adapted to Tiz and T-ix- For any real number s > 0, we denote 

:= {u G H'^(M^) : |x|®^ti G L^{R^)} . 
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According to [12, 5], this Hilbert space can be equipped with the following equivalent 
norms: 




|x|®u||^2 ~ \\'Hl/^u\\l2 + \\n^J^u\\l2 + ||u||^2. 


It is also usefull to recall that, for all 0 < -^ < s, we have 


X 




u 


L 2 


< 




and 




xr’-‘u 


L 2 




( 1 . 6 ) 

(1.7) 


For s > 3/2, is an algebra. Moreover, the self-adjoint operators Tiz and Tix 
generate the groups of isometries 9 i-A and 6 i—>■ on any S®, s > 0. 


To derive the limit model as e —)■ 0, we need to introduce the following nonlinear 
function: 

/ 2cr 

F{e,u) = e^^^^ 


e-m^u 


e-idn^u 


= gi0CH.-l/2) 




2a 




and study the behavior of F (t/e^,u), as e —>■ 0. For s > 3/2, S® is an algebra 
and it is readily seen that T G (7(1^ x S®, S®). Moreover, since the spectrum of the 
quantum harmonic oscillator Hz is {n+ 1/2, n € N}, the operator g*^'(Wz-i/ 2 ) jg 2 ^^ 
periodic with respect to 0, so T is also 27r periodic with respect to 9. Denoting the 
average of this function by 


1 r 

-Fav(tt) := hm 7 ^ F {9,u) 
T —>-00 1 Jq 

2 vryo 


d9 




2a 


,-ienz 


( 1 . 8 ) 


u ) d9, 


the limit model as e —)> 0 reads 
1 . , 1 


idtcj) = --^x4> + 2 (l^l^ “ - ^zLz4> + AFav((/) (1.9) 


with the initial data (j){t = 0) = ipo- The Gross-Pitaevskii type energy associated to 
this equation is 


= 7 ; i \^x4>?‘dxdz + \ [ (|xp - (^23^1 - ^^ 1 X 2 )^) 

2 Jr3 2 7r3 


— il.z{Lz(j), + 
where here and in the following. 


A 


27r(cj -|- 1) 


p A27r 
JR? Jq 


,-ienz 


2 ( 7+2 


4> d9 dx dz, 


{u, v)i 2 = Re 


uvdxdz. 


Note that (1.9) is still a model in three spatial dimensions. Except for the nonlinear 
averaging operator Fav(</>)) however, the variable z only enters as a parameter and 
hence the linear part of the dynamics with respect to 2 : is, in fact, trivial. This 
allows to derive from (1.9) an effective two-dimensional limiting model, provided 
the initial data is polarized on a single mode of Hz, see Corollary 1.2 below. 

One should also note that in (1.9) there is a non-trivial effect due the presence of 
the rotation. Indeed, in the case where the rotation axis is not parallel to the z-axis, 
i.e., if Hi,D 2 7 ^ 0 , a repulsive quadratic potential is present in the limiting model 
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(see also the discussion at the beginning of Section 2.1 below). The reason for this 
effect becomes apparent from the scaling of equation (1.2), which includes a rotation 
term of order The latter becomes large in the limit of strong confinement 

e —>• 0, resulting in an effective centrifugal force counteracting the original trap. In 
the physics literature, it seems that it is almost always assumed that the rotation 
axis is equal to the z-axis, and hence, this effect is almost never considered. We 
finally remark that, at least formally, a second order averaging procedure (similar 
to [4, 10]) can be used to derive (1.9) from (1.2) directly. In order to make this 
procedure rigorous, though, uniform (in e) energy estimates are needed which seem 
to be rather difficult to obtain on the level of (1.2) (given its singular scaling). Thus, 
instead of working with (1.2) directly, we use the change of variables (1.3) which 
yields the same effect as the second order averaging and also allows us to use the 
better behaved model (1.4). 


Our main result is the following theorem. 

Theorem 1.1. Let 1 < u < 2 and ipQ G Then the following holds. 

(i) The limit model (1.9) admits a unique maximal solution (j) G (^([O, Tmax); 5]^) H 
C'^([0,rmax),-f^^), with Tmax G (0,+Oo], SUch that for all t G [0,Tmax)- 

\\4>{t)\\L2 = IIV'oIIlz, E{(j){t)) = Ei'tpo), {'Hz4i{t),(j){t))L2 = (?^^'0o,V'o)l2. 


Moreover, we have the blow-up alternative: 


if Tmax < +00 then lim \\V = +oo. 

i —^^max 


(ii) For all T G ( 0 ,Tmax)) there exists et > 0, Ct > 0 such that, for all e G (Oje-rj) 
(1.2) admits a unique solution G (^([O, T], S^) nC'^([0, T], L^), which is uniformly 
bounded with respect to e & (OjEt] -^°°((0 , T), S^) and satisfies the error bound 


max 

46[0,T] 


fj%t) - 


< Ct £. 

1,2 


As an immediate corollary we have: 

Corollary 1.2. Denote by {Xn, ^n)neN the n-th eigenfunction/eigenvalue-pair of 
the one-dimensional harmonic oscillator T-Lz- Assume that ipQ G T? is such that 

fio{x,z) = ipo{x)Xn{z). 

Then for all T G ( 0 ,Tmax) we have 


max 

iG[0,r] 






T{t)Xr 


L 2 


< Ct e, 


where ip(t, x) solves the effective two-dimensional model 

idtif = -^Az:ip + ^ - (^ 23:1 - ^lX2f'j T “ ^zLzT + T 


( 1 . 10 ) 


with (p{t = 0,x) = ipo{xi,X 2 ) and 


Kn ■■= A 




the effective nonlinear coupling constant in the n-th energy band. 
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This result follows from Theorem 1.1 (ii) and the fact that (1.9) preserves the 
initial polarization, i.e., admits solutions of the form 4>(t,x,z) = (p{t, x)xn{z). To 
see this, recall that the eigenfunctions {Xm}L^ form an orthonormal basis of L^(M). 
Using this we can write 

= Y, e-^^^-Xmiz){Xm,f)L^ 

mGN 

and hence (1.8) implies 

1 _ f^TT 

-F’av(V?Xn) = X] / \Xn\‘^''Xn) ■ 

However, Am — G N and thus, this integral is identically zero unless m = n, for 
which it is equal to 2?!. Thus, the whole sum collapses to one term only and we 
consequently obtain that in the case of polarized solutions, (1.9) reduces to (1.10). 
The latter is an effective two-dimensional model describing the degrees of freedom 
in the unconstrained direction. 

The paper is organized as follows: In Section 2 we shall, as a first step, establish 
well-posedness of the Cauchy problem corresponding to both the three dimensional 
NTS (1.2) and the averaged limiting model (1.9). Once this is done, rigorous error 
estimates between the exact and the approximate solution will be established in 
Section 3. In there, we shall also indicate how to obtain an improved error estimate, 
provided 'ijj^ satisfies sufhciently strong regularity assumptions. Finally, in Section 
4 we shall show how to adapt our results to the situation with strong confinement 
in two spatial dimensions, and derive the associated limiting model. 


2. Analysis of the Cauchy problems 

In this section we shall prove local and global well-posedness results for equation 
(1.2), i.e., the original NTS in d = 3 dimensions, and for the formal limiting model 
(1.9). The analysis of the former is relatively standard and follows along the lines of 
[1]. We shall therefore only sketch the main ideas and rather focus on the Cauchy 
problem corresponding to (1.9). 

2.1. The Cauchy problem corresponding to the averaged NLS model. In 

this subsection we shall analyze the Cauchy problem 

idtcl) = ^ - {^2Xi - (f) - + AFav((/>) (2.1) 

with some general initial data (j){t = 0) = (/>o € Recall that cj) depends on the 
space variables x G and 2: G M, but in this problem dispersive effects only occur 
in the x direction (due to the lack of a second order derivative in z). The basic 
existence proof therefore requires several changes from the standard approach. 

To this end, let us hrst derive Strichartz estimates adapted to the situation at 
hand. We recall that, in dimension two, a pair (q, r) is said to be admissible if 
2 < r < 00 and 

1 _ 1 1 

q 2 r' 
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We denote by U{t) = the strongly continuous group of unitary operators gen¬ 
erated by the Hamiltonian 

H = ^ - ^1x2)^) - ^zLz- (2.2) 

This operator can be seen as a special case of the Weyl-quantization of a (real¬ 
valued) second order polynomial H{x,^). It is thus essentially self-adjoint on 
C'^(M^) C = {/ G : \x\f £ L^(M^)}, cf. [13]. In the case with¬ 

out rotation ilz = 0, H is of the form of an anisotropic harmonic oscillator with 
potential 

V{x)= ^ x\+ ^ xl - ^l^ 2 XlX 2 . (2.3) 

Clearly, this potential becomes repulsive if 121,122 > 1. Physically speaking this 
results in a loss of confinement, and thus, the destruction of the condensate. On 
the other hand, by means of Young’s inequality, one easily sees that a sufficient 
condition for confinement, i.e., V{x) —>■ -|-oo as |x| —>■ 00, is 

12^ -M22 < 1. 

In this case H has only pure point spectrum with no finite accumulation points. 
Similarly, in the case with rotation 12^ 7^ 0 the operator H remains confining pro¬ 
vided Qz is sufficiently small (with respect to I2i, 122). This can be seen by rewriting 
(2.2) in the form of a magnetic Schrodinger operator 

H = ^{V + Aix)Y + V{x)-^nl\x^\^ 

where V is as before and A{x) = QzX'^ with x-^ = (x2, —xi). In this form, the effect 
of the rotation term Lz has been split into Coriolis and centrifugal forces. The latter 
is seen to act as a repulsive quadratic potential, counteracting the confinement. 
Depending on the size of f2i, 122)^22 we thus might have de-confinement due to 
the combined effects of the rotation and the strong conhnement. This also has an 
influence on the question of global existence of solutions to the NLS, see the Remark 
2.3 below. 


Lemma 2.1 (Vectorial Strichartz estimates). There exists 5 > 0 such that the 
following properties hold true. 

(i) For any admissible pair {q,r), there exists C such that, for all (f> £ L^(]R^), 




(2.4) 


where stands/or L^((—<5, (5), L^(M^, L^(M))). 

(a) For any admissible pairs {q,r) and {'y,p), there exists C such that, for all f = 
f{t,x,z), 


' {—S,s)n{s<t} 


U{t — s)f{s)ds 


<C1||/|| 


LlLlLl 


L] Lg L2 


(2.5) 


Proof. This lemma relies on the usual Strichartz estimates for the group U (t) acting 
on functions depending only on x and on Minkowski’s inequality. The existence of 
Strichartz estimates for U{t) thereby follows from the results in [13] which can be 
directly applied to the Hamiltonian given in (2.2). The fact that H in general will 
have eigenvalues, implies that dispersive effects will only be present for small \t\ < 6, 
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preventing the existence of global-in-time Strichartz estimates, cf. [7] for a more 
detailed discussion on this issue. 

Let us introduce a Hilbert basis (ep)pgN of Decomposing the function 

(j) G as (j){x,z) = '^p4’p{x)ep{z), one obtains: 


\\um\ 


LlLlLl 


\\u{m\Li 


< 


’pllLr 


LI 


rr/2 


\\um\h 


r/2 


rq/2 


r9/2 


r9/2 


Lt 


< 


'P\\% 


Here, we have used twice the Minkowski’s inequality [11] in the third line (notice 
that q/2 > 1 and r/2 > 1) and the usual Strichartz estimate for each U{t)(pp in the 
fourth line. This proves (2.4). 

Let us prove (2.5): For f{t,x,z) = Yip denoting 


9 j{t,x) = 


U(t- s)fj{s)ds, 


' {—S,5)n{s<t} 

we estimate similarly (we have again q/2 > 1 and r /2 > 1 ) 


f U{t - s)f{s)ds 

2 

1 X!] 

2 


J (—5,(5)n{s<t} 

LlLlLl 


Ll 

,r/2 

E/x 


r<j /2 


< 


||EI5; 

Eii^i 


r/2 


r9/2 


2 

'piiLr 


r9/2 — 




'pWLr: 


LI 


<c||EIII/r 


Ilp '/2 


<cY.\\s, 


Pll r 7 ' rp' 


cElllli/. 


pi IILp'/2||^7'/2 


Lp'AilY^ 


Here, we have used, in the fourth line, the Strichartz inequality for each gj and, in 
the fifty line, the reverse Minkowski’s inequality [11] (note that we have necessarily 
y/2 < 1 and p'/2 < 1). The proof of Lemma 2.1 is complete. □ 

Proposition 2.2. Let 4>o G S^. Then there exists Tmax £ (0,-|-oo] such that (2.1) 
admits a unique maximal solution (j) G C'([0, Tmax)) E^), in the sense that 

if Fmax <+00 then lim ]]Va;i/'^(t)||i 2 =-7oo. 

Moreover, the following conservation laws hold 

\\m\L^ = UoWl^, = E{cpo), \\ni/^m\L^ = iiwE</>oIIl2. 

Furthermore, if 4>o G T?, then (p G C'([0, Tmax)) s2)ncH[o,r^ax),E2). 
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Note that for (^(t) € we have the strong form of the conservation law, i.e., 

Proof. Step 1: uniqueness in Let u G L°°((0, T), S^), u G L°°((0, T), S^) be 
two weak solutions of ( 2 . 1 ), where, without loss of generality, we can assume that 
0 <T < 6. From (2.1), it follows that 

{u - u){t) =iX f U{t- s) (Fav(^t) - Fi,^{u)) ds 
Jo 

and then, by (2.5), for any admissible pairs {q,r) and ( 7 ,p). 


11'*^ < Cl ||F’av(^) Lf L'^ ' 

Let us choose r > 2 and p > 2 such that y + ^ = 1 ~ f) where we recall that 
1 < (T < 2. Then, denoting v = and v = one gets by Holder, 


l-C'av('w) Tav(tt) ll^p'^2 — 2 y]- 


< C 


< c 


< c 


L 

p27r 

/ 11(1 

Jo 


p/ d6 
Lx 


V]2a _ u|||^p/^2 de 


rZTT 

/ + ll^lli4L-)ll^ - ^hrLlde 

^ 0 
p27V 

/ {\\vf^r + M^^r)\\v-^\Lr^dO 
Jo 


= c’(ii«iii'^i+Hum'll) iiu - 


Here, we have used the unitarity of in and and the embeddings 

iCi(M3) L^(M 2 ^Lf'(M)) and ^ Ll(R‘^, Ll{R)) (see the Appendix of 

[ 6 ]). Since u and u belong to L°°((0, T), S^), this yields 


\U — U\\iq]^rl2 < C\\U — U| 




Since 7 ' < 2 < < 7 , this inequality is enough to conclude that u = u, see Lemma 4.2.2 
in [9]. 


Step 2: local existence. Let us adapt to (2.1) the proof of well-posedness of NLS 
with a quadratic potential of [7], using Kato’s strategy (see for instance [9]) and the 
vectorial Strichartz estimates given in Lemma 2.1. The main technical difficulty 
here is the fact that we are working with an NLS in three spatial dimensions, but 
we can only utilize the dispersive properties of the two-dimensional Schrodinger 
group U{t) = In order to remedy this, an important ingredient will be the 

anisotropic Sobolev imbeddings proved in [ 6 ]. 

With this in mind, local in-time existence for solutions in (/)(t) G can be proved 
by means of a fixed point theorem applied to Duhamel’s representation of (1.9), i.e., 

f){t) = U{t)(j)o -iX [ U{t- s)Fs,v{(l)){s)ds =: 4>((?i>)(t), 

Jo 
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where, as before, U{t) = and H is given by (2.2). We want to show that for 
4>0 ^ and sufficiently small T > 0, $ is a contraction mapping in the complete 
metric space 

IIV'lli,j°°(si) + IIV'IIl«lji, 2 + + l|vV'lli,5i,jL2 < m}, 

equipped with the distance 

d{u,v) = \\u-v\\l^l 2 l 2 + \\u-v\\^LiL 2 . 

The real numbers r, q, p, 7 are taken as in Step 1 above and T, M are to be chosen 
later. To prove that Xt^m is stable by <f>, one first checks that the commutator 

[d,,H] = [z,H]=0, 

whereas 

[x, H'\ — XX ; 

with x~^ = (xi, —X 2 )- Similarly, we find 

lVx,ff] = VxV + ii 2 zV^ 


where V is as in (2.3), and hence XxV is in fact linear in x. We consequently obtain 
that the combination of (f), x(f) and X(j) form a closed coupled system of equations. 
Therefore, by applying the operators x and V to (2.1) and by using Lemma 2.1, we 
obtain for all (p G Xt,m, 


^ II<('o||l 2 + ||x0o||i2 + ||V())o||l 2 + ||x(/)||^i^2^2 + \\XxP\\lIlIli 

+ lp' 1,2 + 

< ||<).o||si+TM+||Fav((/.)|| 


L7'Lg'l,2 + ll^^av(</')||^7'^g'^2 ; 


il7'l£'l2 


+ ||xFav((/>)|| 


l 7 'l£'l 2 


+ l|VFav(</.)|| 


l 7 'l£'l 2 


where a <b stands for a < Cb for some constant C > 0. Denoting v = one 

gets 

r 2 iT 


I V.Fav (</>)!!, 


< 


< 

r\j 


< 


1 

/*27r 

/ \\v\\%L-°\\^M\L-iLldO 

Jo 

p27V 

/ ll^ll|"l|V.u||i.i2d0 = 27r||0|||^||V, 
Jo 


where we have used the fact that is unitary in ^^^1 together with a 

Holder estimate and the embedding Lf^(M^) L^(M^, L“(M)). Hence, 


I V^Fav (</>)!! 




Q-l' 


< TXT 




Similarly, we obtain 
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To estimate zFi^^{u) and dzFa;y{u), we use (1.6) several times: 


\\^Pa.v{'P)\\j^p' 1^2 

+ \\dzFav{p)^^i^ 








r2'K 



/*27r 



< 

rs_/ 


1 (|u 

\^^v) de 

= 

/ 

(|u 

de 




0 

LiLl 

Jo 


/*27r 




/*27r 




^ J 

( 


, de 

iS LI 

s/ II 

Jo 

\v\‘^‘^{\zv\ 

+ l-^wl) 

l£'l 2 de 


r2'K 








~ J 

( 

V 

\%Lf> 

LMl + II 

dzVhrLl 

) de 




/*27r 








^ J 

( 

V 


x^LldO = 



r 1^2 
x-^z 





which yields, again 

Finally, we have proved that 


l) ’ 


y' T p' 
-'t 




_ / 

< C\\(j)ohi + CTM + CT^ 

We now set 

M = 2C||</.o||si 

and choose T small enough so that 

Q-l' 9 ,1 M 
CTM + . 

It follows that $((/>) G Xt^m- The contraction property can then be proved by 
following the same lines as the proof of uniqueness in Step 1: it can be proved that, 
for T small enough, we have 

for all £ Xt,m- Hence, by Banach’s fixed point theorem, <h has a unique fixed 
point, which is a mild solution of (2.1). 

Step 3: blow-up alternative. From the uniqueness result and from the fact that 
the existence time in Step 1 only depends on ||i?io||si) define the maximal 

solution (p G (^([O, Tmax); 5]^) and obtain a first blow-up alternative in terms of the 
whole norm: 


if Tmax < +00 then 


lim 


.1 = + 00 . 


Then, we compute from (2.1) 


d 

dt 


|3;</’(^)lli2 = 2Im / X ■Vx<p{'t,x) 4>{t,x) dxdz < \\x(l){t)\\l 2 + \\Xx<p{'t)\\‘i 2 . 


Hence, a bound on ||Va;(/>||L2 yields a bound on ||x(/)(t)||i2 by the Gronwall lemma. 
Since the norm of (p is conserved, it is clear that limt_>.7’^^^ ll'/’(f)llsi = +oo 














DIMENSION REDUCTION FOR ROTATING EEC 


11 


implies that ||V0(t)||i2 = +oo. We have proved the blow-up alternative 

as it is stated in the Proposition. 


Step 4- conservation laws. In order to prove the conservation laws stated above, it 
is enough to consider the case of local-in-time solutions (f){t) which are sufficiently 
smooth and decaying. By following a standard regularization procedure (as given 
in, e.g., [9]), this can then be extended to general (j){t) G Conservation of mass 
then follows from the fact that H is self-adjoint and hence 

^II<('WIIl 2 =Re AFav((/>),0))L2^ = 2AIm(Fav(0),(?i))L2. 


However, 


Im(Fav(((')»)L2 = 







L2 


1 






de = 0 , 


which implies \\(t){t )\\\2 = II<(' 0 II12- 

A similar argument can be used to prove that ||?^y^(^(t)||^2 = ||?^y^(()o||i,2, having 
in mind that both Tiz and H are self-adjoint, and that the following commutation 
relations hold: = 0 as well as = 0, see [5]. 

Finally, in order to prove the conservation of the energy it is useful to first note 
that E{4>) is formally equal to 


E{(l)) = + X [ Ns;^{(j))dxdz, 

Jr3 


where 


2^(a + 1)1 


2n 


..-idHz 


2a+2 


de. 


We then compute 
d 


— = 2 Im(Fav(</>),-R</>) l 2 = 2Im(Fav(((>),*5t(/> - AFav(<(>))L2 

= -2ARe(Fav(^),5t(())x,2 


From here, it follows that E{4>) is conserved provided that 

1 p 

2Re(Fav(0),9t(())x,2 = — / N^^{(l))dxdz. 

dt Jr3 

This, however, can be seen by a direct computation, using the definition of T'av 
-^av? i-G., 


Re(Fav(</>),5t(())^2 




de 



2a 

dt 


e-^ 0 H. 


2 

de dx dz 


1 d 
47r(cj -I- 1) dt 





2a+2 

de dx dz. 


Step 5: E? regularity. Assume that (f)Q G Since is an algebra and 
is unitary on it is easy to see that Rav is locally Lipschitz continuous on 
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provided <7 > 1. Hence, by a standard fixed-point technique, we can show the 
existence of a unique maximal solution (p G C([0,Ti),S^) H (^^([O,Ti),L^), with 
0 < Ti < Let us prove that Ti = by contradiction. To this end, we 

assume that Ti < T^ax and consequently deduce that 


lim ||</>(t)||s 2 =+00 and sup ||(?!>(t)||si <+oo. 

tG[ 0 ,Ti] 

Therefore, it suffices to find r < Ti such that ||<^(t)||s 2 is bounded on the interval 
(r, Ti) to have the desired contradiction. Let r > 0, to be fixed later, be such that 
0 < Ti — T < (5 (5 is defined in Lemma 2.1). By differentiating (2.1) with respect to 
time, we obtain that 

dtpit) = U{t)dt(p{T) J U{t- s)dtFs_y{(j)){s)ds. 

Denoting pj = we compute now 

dtF^AP) = ^ dti?) d6 + ^ dO 

and, using the same admissible pairs (g, r) and ( 7 ,/?) as in Step 1, we estimate 
similarly 

Using the Strichartz estimates (2.4) and (2.5), this yields 


\\dt(t>\\LlL^^Ll < C!\\dt<p{T)\\L2 +C\\dt(, 




where the time integral is computed on the interval I = (Ti — r, Ti). Since 7 ^ < g, 
if Ti — T is small enough, this implies 


\\dt4>\\LlLr^Ll < C\\dt(p{T)\\L2 
and, applying again the Strichartz estimates, we get 


\\dt(p{t)\\L'^Ll^^ < C'||(9t(/>(r)||i2 < C'||(/>(t)||s2. 

For the last inequality, we used (2.1) at t = r and the following estimate deduced 
from a Sobolev inequality and from interpolation inequality 

p27V 

WFMWl^ < C / < C\\u\\l-+^ < C||u||“.||u||^"+'-“, ( 2 . 6 ) 

Jo 

with s = 2 a+i ^ ~ max ((7 — 1,0) < 1. 

Now, we remark that, by integrations by parts, a direct calculation gives 

WLzpf = -{xlp^dl^p) - {xlp,dl^(p) + 2Re{xiX2p,d:cida,^4>) -^(^^12 
so, using directly ( 2 . 1 ), we estimate 

||Aa;(^||x,2 < 2||5t0||j;^2 -b C'|||xp0||j;^2 + C\\Lz(l)\\i,2 + ||Fav (0) ||i;,2 
<2\\dML^+C\\\x\^P\\L2+C\\/2.z,P\\^ff\\\x\^^^^ 
and then, using the bounds of ||(/>||si and || 5 t(/>||^ 2 , we deduce that, for all t G / 

II^3;<?!>(^)||l 2 < C-b C'||(/>(r) 11^2 + C'|||xp(/>(t) 11^2. (2.7) 
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Next we can proceed similarly as above to estimate ||| 3 : We have 

\x\^ct){t) = U{t){\x\^cj){T)) + A J U{t - s) {2Ax(p + 2x ■ SIx<p + |xpFav((?!))) {s)ds. 
Hence, using that 

and that, by (2.7) and (1.7), 

IIAj;^ + X ■ Va;(/>||i 2 <C + C'||(/)(t)||s 2 + C'|||3:|^<(>||l 2L2, 
we get by Strichartz estimates 

llkP<('llLf=LiLi + llkP<('llLfLJi2 ^ C'^H-||(()(r)||s2 + |||xp0||j;^ij;^2x,2 + |||x|^0||^y^^^2^ 
From 7 ' < q, it is easy to conclude that, for Ti — r small enough, 

\\\A^4>\\l'^lIlI + llkP</’lli,|Lji ,2 <C + C'||(()(r)||s 2 
and, with (2.7), that 


W^x’PWbf'Ll^^ <C + C'||((>(r)||s 2 . 

Finally, consider the equation satisfied by T-iz4>- We have 

Uzckit) = U{t){'Hz(t){T)) + \ J U{t - s)'HzF^^{(t)){s)ds 

and, denoting again ■0 = one computes 

r 2 -K 


( 2 . 8 ) 

(2.9) 


||H.Fav(0)||L2 <C 

r 2 -K 


0 


+ k|'| 0 |'"+'|L 2 dO 


< c 


'0 


(V’IIli + wnfF^wdMi^ + wmuwzmLi do 


< c 


1 ‘ 2 'k 

Jo 


2(7 


ILJ° WF-z'^PWL'^dO, 

where we used the following Gagliardo-Nirenberg inequality in dimension 1: 

115^1x11^4 < C'||u||^2^||u||]/i < C\\nzu\\]^2\\u\\L^- 
Hence, as above, we get 


P.i"av(0)||^7^2<C' 


1*277 

Jo 


"^iWnMLrLidd = cumnzPhrLi 


2(7 I 


which enables to conclude again with Strichartz inequalities that, for Ti — r small 
enough, 

WTdzPWl^LlLl + \\Tdz4>\\L^LlLl < C + C\\(I){t)\\y.2 . ( 2 . 10 ) 

From (2.8), (2.9) and (2.10), we deduce that ||0(t)||s2 is uniformly bounded on the 
interval I = (r, Ti): the proof is complete. □ 


Remark 2.3. There are certainly situations for which Tmax = + 00 . In view of 
the discussion at the beginning of Section 2.1, this will be true, in particular, if 
nf + nl < 1 and the nonlinearity is defocusing A > 0, since in this case, the results 
given in [1] apply. However, if the effective centrifugal force is too big, the resulting 
repulsive quadratic potential requires particular techniques, see [ 8 ] which would need 
to be combined with the effect of the rotation term. 
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2.2. The Cauchy problem for the NLS equation in 3D. Before we can proceed 
to the proof of convergence for solutions of (1.2) as e —>■ 0, we need, as a 
final preparatory step, a suitable existence result for the three-dimensional Cauchy 
problem (1.2) 

Proposition 2.4. Let ipQ ^ Tj^ . Then, for any fixed e > 0, there exists Tf G (0, -|-oo] 
such that (1.2) admits a unique maximal solution G (^([O, Tf), S^), i.e., 

if Tf < -|-oo then lim^ I|VV'^(^)IIl2 = +oo. 

Furthermore, if'ifo G 5]*, s > 1, then G (^([O, Tf), S^) n (^^([O, Tf), L^). 

Proof. The proof follows along the same lines as the one for Proposition 2.2 above, 
i.e., through a fixed point argument for Duhamel’s formula in a suitable metric space. 
Indeed, it is even easier, since (1.2) is a standard three-dimensional NLS equation 
with quadratic potential and rotation term. For such equations, the local and global 
existence theory in based on Strichartz estimates and the use of energy-met hods, 
has been studied in detail in [1] (see also [7]). Additional smoothness for ■00 £ 51* 
with s > 1, then follows by the same arguments as given in Step 5 in the proof of 
Proposition 2.2. □ 

One might be concerned that, as e —>■ 0, the existence time Tf 0, but our 
convergence proof below will show that this is indeed not the case. We finally, note 
that (1.2) admits the usual conservation laws for the mass and the total energy. 
The latter, however, is in general indefinite, due to the appearance of the angular 
momentum operator. Since, we shall not use any of these conservation laws in the 
following, we omit a more detailed discussion of these issues and refer the reader to 
[!]• 


3. Convergence proof and error estimates 

This section is devoted to the proof of our main Theorem 1.1. Item (i) of this 
theorem is a consequence of Proposition 2.2.We prove items (ii) in Subsection 3.2 
after we have obtained some uniform estimates. 

3.1. Uniform estimates. Let 0 < T < Tmax- By Proposition 2.2, we know that 
the solution (f of (1.9) belongs to C'([0,T],S^) H C'^([0,T], L^). Using the Sobolev 
embedding > L°°(M^) and the unitarity of in we have 

||g itnz/e 0||^oo((o_7’)xr 3) < Clle 0 ||l°°((o,t),e 2 ) = C'||^||x,oo((o_t),e 2 ) < +oo, 

so the following quantity is finite; 

M := sup ||e"*‘'^"/^VllL-((o,r)xK3). (3.1) 

£>0 

Notice that, in particular, we have HV'oll LjOO — ll<^(t = o)|| ]^oo < M. 

By Proposition 2.4, (1.2) admits a unique maximal solution ■0'^ G (^([0, Tf), S^) n 
C'^([0,Tf),L^). We recall that the function defined by 

yf — g-i£z{QiX2-n2Xi).^e 

satishes (1.4). By (1.7), it is clear that, for all t G [0,Tf), 

(1 - C'ie)||V'^(t)||s2 < ||^x^(t)||s2 < ( 1 + C2e)||i/’^(t)||s2- 
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Moreover, since 


.(fiix2-n2xi) _ ^ -i£z{niX2 - n2Xi) r 

Jo 


ds, 


we have also that 

Hit)-^l>%t)U.<Ce\\u%t)\\^.. (3.2) 

This will allow us to infer the desired approximation result for once we have a 
sufficiently good estimate on the difference between and the limit (j). 

From a Gagliardo-Nirenberg inequality, we get 

\\ul - V’oIIl- < C\\ul - - Mh2 < 

Hence, for e <e\\= (M/2C'i)^, we have 

||ii^(0 )||l°° < ||iio “ V'o||i,°° + ||^o||l°° < 3M/2 

and we can define 

= sup{t € [0,Tf) : for all s G [0,t], ||u^(s)||loo < 2M} . (3.3) 

Lemma 3.1. There exists a constant Cm such that, for 0 < e < ei and for all 
t € [0, T^], we have 

||ii^(t)||s2 < Cm- 

Proof. We first recall from (1.6), that 

Il^lls2 — II^2 '^IIl2 + ll^a;^llL2 + ||ll||/^2. 

We will now derive suitable a priori estimates for these three parts of the S^-norm. 
To this end, we hrst multiply (1.4) by u^, integrate over R^, and take the real part 
of the resulting expression. This yields 


d 11 c-119 


i.e., the conservation of mass. For the other two parts of the S^-norm, we hrst 
compute the commutation relations 

[Hz, TLx] = [Hz, Lz] = [Hx, Lz] = [Hz, fllXi ± 1123^2] = 0, 

(where in the third expression we have used that Tix is rotationally symmetric), as 
well as 

[Hz,z] = -dz, and [Hz,z^] = -{l + 2zdz). 

Keeping these relations in mind, we can thus apply Tiz to (1.4), commute, and, 
after multiplying by TizU, integrate over R^. Taking the real part of the resulting 
expression, we obtain 
d 


■^\\H.zu'^\\l2 = + Q2 ) 1'<m{'HzU^, + ‘^zdzU^'iii 

+ 2 elm{'HzU^,Tlz{^iXi + Q.2X2)dzU^ - 2{Q.2dl^^^u^ - ^2 


+ 2Xlm{'HzU^ ,Tiz\u^\^'' u^) L'^ ■ 

After several integrations by parts and using Cauchy-Schwarz, this yields the fol¬ 
lowing estimate 

^ \[ Hzu ''\\\2 < e^Ci {\[ Hzu ''\\\2 + \\ u ^\\\ 2 ) + eC 2 + \[ Hxu ''\\\ 2 ) 

+ 11^2 + , 
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where Ci, C 2 , C3 are some constants depending only on fii, O2, and but not on 
e. Similarly, a computation shows 

■^W'HxU^\\'j ^2 = ^^2 + “^xidx-^u^)lni{'HxU^, + 2 x 2 dx 2 U^) 

— 2 ^ 71^2 lina{'HxU^, X2dx-^v!^ + xidx2 U^)l'^ 

+ 2 eQ.z lm.{'HxU^, z{VlidxiU^ + ^29x2^^ — (^^12:2 — ^ 2 Xi)vF)) 

+ 2 \lm.{T-LxU^ ,T-Lx\u^^'^u^)l^, 

and using again Cauchy-Schwarz yields the analogous estimate for ||?^2;U^||^2. Com¬ 
bining the three estimates obtained above, allows us to write 

d 


dt 


u 


Is 2 < Ki \\ u '^\\^2 + \X\K2\\\u'^\ um ^2 


(3.4) 


where i^ri,2 = Ki^ 2 {£, XI 2 , Xlz) > 0 are both bounded as e —)■ 0. Now, we use the 
fact that, by Sobolev’s imbedding, 

equation (3.4) implies 
d 


dt 


\u^\i ^2 < Ki\\u ^\\^^2 + 


IS2- 


Gronwall’s lemma consequently implies that ||ir^||s2 stays bounded for all t € [0, T^]. 

□ 

We note that an important consequence of this lemma is that 

if r*" <-hoo thenT^<Tf and ||u^(T'')||loo = 2M. (3.5) 


3.2. Proof of the error estimate. In this section, we prove Item (ii) of Theorem 
1.1. Consider the function u^. This function satisfies the equation 

idtv^ = -^Axv‘^ + ^ - {Q2X1 -^1X2)^^ -QzLzV^^+XF -t-eri-t-e^r| 

with 

v^t = 0) = 

and where we have denoted 

rf = J, J := {2{VL2dxi — Xlidx 2 ) — Xlz (flixi -|- XI 2 X 2 )), 

We have G (^([O, Tf], S^) fl (^^([O, Tf], L^) and, by Lemma 3.1, 

max ||u^(t)||s2 = max l|ri^(t)||s2 < Cm- 
iG[o,T®] te[o,r®] 

Hence, by (1.7), we get for t G [0,T^] 

\l 2 < 

,£|l ^ _ 


< = C||.2e- —“ v^w^i 


^Clle- 


^^Ils2 = Clll^^lls^ < Cm, 


(3.6) 
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< = Cm\j,2 < Cm. (3.7) 

We are now ready to estimate the difference w'^{t) = v'^{t) — (j)(t) for 0 < t < 
min(r, T'"). This function satisfies 

W^{t) = U{t){ul - V'o) + A U{t - s) (^F - F ds 

+ X J U{t- s) (^F (^, </>('§)) - Fav((/>(s))^ ds 

+ ef U{t - s){rl{s) + £r 2 {s))ds 

Jo 

= + ^2 + ^3 + A 4 . 

Since U(t) is unitary on L^, (3.2) yields 


Pi||l 2 = ||u^ - V’o||l2 < Ce. 

Moreover, for 0 < t < min(T, T*^), (3.1) and (3.3) imply that 


i2||L2 <C [' |u"(s)|2'^u"(s) - |e-“^^/"V(s)P"e-“^^/"V(s) 
Jo 


ds 


< C 


S ) 11 ‘j^ca + 




and (3.6) and (3.7) give 


I (“f 

/ I|w^(s)IIl2c?s 
Jo 




(j){s) 


2a 


i|u^(s) — e 


L2 

—isHz/e^ 


(j){s)\\L 2 ds 


I|714||l 2 < Ce. 

Let us estimate ^ 3 . To this aim, we introduce the following function, defined on 
M X S2, 

F{9,u)=[ {F{s,u) - Fi,^{u))ds. 

Jo 

Since F{.,u) is 27r-periodic and Tav is its average, F{9,u) is periodic with respect 
to 9. ffence, it is readily seen that this function satishes the following properties: 

if ||tt||s 2 < i? then sup 117^(0, u)||s 2 < 

0GR 

if ||tt||s 2 + \\v\\l'^ < R then sup ||i7uJ'(0, tt) ■ v \\]^2 < 

0gR 

Recall that U{t) = where H is the Hamiltonian defined by (2.2). Hence 


U{t-s) (^F (^,0(s)) - Tav((/>(s))) 

= (u{t - s)F (^)‘?^('®))) + - s)HF 

- e‘^U{t - s)DuF -dtcpis), 
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and then 


P3||L^<e"|A| 

+ e^\, 

< Cs^, 


T 


>(i)) 


L 2 


ds 


\ f \\DuJ^ -dMs] 


L 2 


ds 


where we used that (j) G L°°([0, T], S^) and dt4> G L°°{[0,T], L'^). In summary, we 
have proved that, for all t G [0, min(r, Tg)], 


iv^{t)\\i 2 <C£ + C 



12 ds. 


Thus, Gronwall’s lemma yields 

= ||^"(t) - m\L^ = lk"(i)llL2 < Ce. (3.8) 

In particular, we deduce from (3.1), from a Gagliardo-Nirenberg inequality, from 
(3.8) and from Lemma 3.1 that 


u^{t)\\L^ < M + ||u'^(t) - e *‘^^/^V(i)l|L°° 

< M + ||u^(t) - e-*‘«^/^V(i)llii"l|u"(t) - 

<M + Ce^/^i\\u%t)\\j,2 + \\mk^f" 

< M + Ce^/^ 


Hence, for e < et '■= (M/2C')^, we have 

Vt < min(T, T^), ||u^(t)||ioo < 3M/2. (3.9) 

It is clear then that > T. Indeed, otherwise this would imply that < +oo 
thus, by (3.5), that ||u'^(T£)|| = 2M, which contradicts (3.9). Consequently, (3.8) is 
valid on [0,T] and, together with (3.2) and Lemma 3.1, this yields 

_ e-*‘«^/^V(t)llL2 < writ) - u%t)\\L2 + ||u^(t) - e-*‘«^/^V(i)llL2 < Ce 

for all t G [0,T]. We have proved Item (ii) of Theorem 1.1. □ 


We note that under under sufficiently high regularity assumptions on rit), a 
slightly stronger approximation result can be proved. In this case, one can show 
that, indeed, ||H 4 ||j ;^2 < Ce‘^, and not only 0(e) as shown above. To see this, we 
expand 

Rl := £ f U{t — s)rl{s) ds, 

Jo 

using the eigenfunctions of T-Lz- Writing v^{t, x, z) = x)Xmiz) we obtain 

Rl = jf u{t-s) 

mT^neN 

where J is as above. Here, we have also used that is odd and thus only 

indices m ^ n appear in the double sum above and Rf is consequently seen to 
be highly oscillatory. After an integration by parts in time (which requires the 
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improved regularity of one obtains that R\ = Using this one can show 

that, for 0 < T < T^ax and 0 < £ < et, the following improved error estimate holds 


max 

iG[0,T] 




L2 


<CTe\ 


where (j)^ is the solution of (1.9) with initial data 

= 0 ) = 


(3.10) 


Note, however, that if we apply this e-correction to the Cauchy data, the solution 
(p^ does not remain polarized any more, in which case the reduced model is still 
posed in three spatial dimensions. 


4. The case of strong two-dimensional confinement 


In this section, we briefly discuss how to obtain the limiting model in the case of 
strong two-dimensional confinement within the original (three-dimensional) Gross- 
Pitaevskii equation. To this end, we we start with the analog of (1.1), given by 

idtip = -^^'4’ + iU • (x A V) V' -b (4.1) 

subject to 'ip{t = 0 ,x) = £~^'iPq{x/£, z), where as before x = (x,z) € with 
X = {xi,X 2 ) G and z G M. Note that in comparison to (1.1) the roles of x 
and z have been reversed. Thus, in (4.1), the x variables are now the ones which 
represent the strongly confined directions, and we consequently aim to derive an 
effective model depending on the z-variable only. To this end, we rescale 

x'=—, z'= z, 'ip‘^{t,x', z') = eip (t,£x', z') 

and assume that i.e., an even weaker interaction regime as before. The 

rescaled NLS equation then becomes 

idt'ip^ = - ^-z - VL2dxi'>P^) - 

— ie (02X1 — O 1 X 2 ) dzip^ + (4.2) 


with ?/i^(t = 0,x,z) = iI)q{x,z). In order to get rid of the singular perturbation 
we invoke the same change of variables (up to a sign) as in the case of a strong 
one-directional confinement, i.e., 

iA"(t, X, z) = 

After a somewhat lengthy computation, this yields the following analog of (1.4): 

idtu^^ = ^TixU^ + R-zu’^ — -(f^i + — QzLzU^ + 

^ 3 (4.3) 

-b eVLz (Oixi -b O2X2) zu’^ -b 2ze(OiX2 — ^2Xi)dzu‘^ + -e^(02Xi — OiX2)^u^. 


In order to average out the fast oscillations stemming from we introduce 

2(t 


G(0,u) = e 


iOH:, 


e 


e 


which satisfies G G C(R x S®;S®). Moreover, G is easily seen to be a 27r-periodic 
function in 9, since the spectrum of the two-dimensional harmonic oscillator is 
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given by {A„ = n + 1 , n G Nq}- Note however, that the eigenspace corresponding 
to Xn is (n + l)-fold degenerate. We consequently denote the associated averaged 
nonlinearity by 



and find, that, as e ^ 0 , the new limiting model becomes 


de, 


idt4> — ~2 2 ~ ^ 2 )) ~ ^zLzfp + ^Giiv{(j)) ( 4 . 4 ) 


Again, we note the appearance of an additional negative (repulsive) quadratic po¬ 
tential, provided ^ 1,02 7 ^ 0 . 


The limiting model (4.4) has the drawback to still be an equation in three di¬ 
mensions. But, having in mind that [Hx,Lz\ = 0, there exists a joint orthonormal 
basis of eigenfunctions {XalaeN^ where a = ( 01 , 02 )) such that 

^zXct — tJ'aXct: TdxXa — ^nXa: G 1^) ^ — Oil T ®2- 

The simplest situation is then obtained for initial data (f>Q concentrated in the 
eigenspace corresponding to the ground state energy Aq = 1. This eigenvalue 
is known to be non-degenerate, i.e., (j)o{xi, X 2 , z) = ipo{z)xo{xi, X 2 )- In addition, 
Xo = Xo,o is known to be radially symmetric which implies /ig = 0. By the same 
arguments as earlier (see the remarks below Corollary 1.2), we consequently obtain 
that (4.4) admits polarized solutions of the form 

(p{t, X, z) = z)xoixi,X 2 ), 

where (p solves the one-dimensional NLS equation 

dlp> + ^{l- {nl + z^(p + (4.5) 

where kq = A||xo|li 2 ai- 


Our main result in this section is then as follows: 


Theorem 4.1. Let 1 < a < 2 and G Then the following holds: 

(i) The limit model (4.4) admits a unique maximal solution cf G (^([O, Tmax)) S^) O 
C^{[0,Tjnax),L‘^), with Tmax G (0,-t-Oo], SUch that for all t G [0,rmax)- 

\\4>{t)\\L^ = WMl^, E{(f){t)) = E{lpo), {'Hx(j)it),(p{t))L2 = {'HxfJo,ifo)L^. 


(ii) For all T G (0,Tmax)) there exists et > 0, Ct > 0 such that, for all e G (OjEt]) 
(4.2) admits a unique solution ip'^ G C([0, T], S^) nC^([0, T], L^), which is uniformly 
bounded with respect to e ^ (0,eT] in T°°{{Q,T),T?) and satisfies the error bound 


max 

iG[ 0 ,T] 




L 2 


< Ct £■ 


(hi) If moreover the initial data is such that fiQ{x,z) = (Pq{z)xo{x), then, for all 
T G (0,Tmax), wo have 


max 

tg[ 0 ,T] 


fi^{t)-e-^^/^\{t)xa 


< Ct£, 


where ip{t,x) solves (4.5). 
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Proof. The approximation proof follows along the same lines as the one for Theorem 
1.1, up to adjusting the notation (i.e., switching the roles of x and z). The main 
difference concerns the proof of well-posedness for the limiting equation (4.4). In 
contrast to the case of one-dimensional confinement, equation (4.4) only admits 
dispersive properties only in one direction, which might not be sufficient for the use 
of Strichartz estimates. However, since we are working in ^ L°°(M^,j), local 
in-time well-posedness follows from standard arguments, see [9]. □ 

In comparison to initial data polarized along the ground state Aq, the situation 
for initial data polarized along some higher energy eigenvalue Xn, n > 1, is much 
more complicated, due to their (n-|-l)-fold degeneracy. The corresponding solutions 
are then of the form 

(/>(t, X, z) = ^ 

ai+a2=n 

where the coefficients ip a = Pai,a 2 solve a system of n -|- 1 coupled NLS. The latter 
mixes the (pa through the nonlinearity and describes the dynamics within the n-th 
eigenspace. The precise form of the NLS system is rather complicated and hence, 
we shall leave its details to the reader, in particular, since one anyway might prefer 
the description of the dynamics via the effective model (4.4) instead. 
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